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ROTABLE INVENTORY CALCULATION METHOD 

TECHNICAL FIELD 

The present invention relates to a method for conducting repairable inventory analysis, 
and more particularly to a method for computing repairable inventory requirements under 
uncertainties. 

BACKGROUND OF THE INVENTION 

Many industries use repairable, or "rotable," inventories for economic reasons. 
Rotable parts are different from "expendable" parts, which are parts having a low enough 
value that the repair of such parts does not make economic sense. Rather, such expendible 
parts are merely discarded and replaced with new parts. Rotable parts, by contrast, tend to be 
more expensive, making their repair and reinstallation, rather than simple replacement with a 
new part, more economically justifiable. 

Repairing a part or getting a new replacement can have uncertain lead times (i.e., the 
time interval between ordering a replacement and its delivery), and therefore rotable part 
inventories are used to bridge the gap between demand for the part and its supply as well as to 
maintain a selected high level of customer service. Due to the high cost of rotable parts, 
however, it is desirable to minimize the number of rotable parts held in inventory. But 
balancing minimal inventory with a desired customer service level (i.e., a measure of 
customer service defined as a ratio of parts delivered on time to the number of parts ordered) 
is difficult because the lead times for repairing parts and obtaining replacements are uncertain. 
Thus, there is a desire to calculate a minimum rotable inventory level that can satisfy a given 
customer service level. 

Currently-known attempts to create models solving this problem have not provided 
satisfactory solutions because they use deterministic methods that assume parts arrive into a 
repair shop and are repaired according to standard time distributions. In actual practice, 
however, arrivals and repair time are much more uncertain. Thus, current models for 
optimizing rotable inventory do not generate satisfactory solutions because they fail to take 
these uncertainties into account. Further, assets may contain multiple types of rotable parts 
and the asset service level is determined by the service levels of individual part types. The 
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interrelationships between those part types make it more difficult to determine optimal 
inventory levels for individual part types to achieve a desired assert service level. 
Deterministic methods are unable to consider the interrelationships between the parts. 
[5] There is a desire for a method that can calculate the optimum amount of rotable 

inventory needed to satisfy a given customer service level while taking these uncertainties 
into account. 

SUMMARY OF THE INVENTION 

[6] The present invention is directed to a method for conducting a repairable inventory 

analysis that calculates the minimum number of repairable parts to be held in inventory while 
still being able to satisfy the demand for spare parts even in view of uncertain repair lead 
times. Given a service level, which is defined as a total number of on-time delivered parts 
delivered by the total number of requested parts, one embodiment of the method uses a set of 
equations that relates the optimal rotable inventory level to a customer service level, given 
arbitrary probability distributions of a repair lead-time and a part arrival process. A search 
procedure is used to find the optimal rotable inventory level. 

[7] In one embodiment, it is assumed that the number of parts requiring repair in a given 

time unit is constant and a part repair lead-time distribution is given. The mean and variance 
of the number of available parts at any time (or the amount of back orders, if the mean is 
negative) are computed. In this embodiment, this number is random and cab be proven to be 
normally distributed. A probability term that corresponds to a given probability level under a 
normal distribution is computed or looked up from a standard normal distribution table. The 
corresponding inventory level is computed from the probability term, the mean, and the 
variance. The number of on-time deliveries is then calculated based on its distribution, which 
is obtained by truncating the above-mentioned normal distribution with [0,A] and lumping the 
probability in each tail to the corresponding boundary. If the on-time delivery value is 
acceptably close to a desired on-time delivery value, the inventory level is deemed to be 
optimal and can meet desired customer service levels without having excessive, costly extra 
parts on hand. Otherwise, the probability level is updated following a search procedure, and 
the above calculation is repeated until the optimal inventory level is obtained. 
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[8] In another embodiment, it is assumed that the number of arriving parts in a given time 

unit is random. The algorithm described in the previous embodiment still applies, except that 
the distribution of on-time delivery is more complicated and difficult to obtain in a closed 
form. Equations for computing the above distribution numerically, and for computing the 
mean and variance of the number of available parts at any time, are provided. 

[9] In another embodiment, an asset with multiple part types is analyzed. The asset service 

level is computed as the product of service levels of individual part types. With the above 
algorithm, inventory levels for individual part types can be obtained with given service levels, 
and the problem is formulated as finding the optimal service levels for individual part types. 
With known per-unit inventory costs, the optimization problem is formulated as minimizing 
total inventory cost, subject to the above service level relations, and the inventory-service 
level relationships for individual part types. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[10] Figure 1 is a block diagram illustrating a rotable inventory problem to be addressed by 

the invention; 

[11] Figure 2 is a graph depicting examples of on-time delivery distributions; 

[12] Figure 3 is a flow diagram illustrating an example of a method used to find an optimal 

inventory level according to one embodiment of the invention; and 
[13] Figure 4 is a block diagram illustrating one embodiment of a computer system that can 

implement the inventive method. 

DETAILED DESCRIPTION OF THE EMBODIMENTS 

[14] The invention determines a minimum rotable inventory level that can still satisfy a 

given desired level of customer service by generating a set of equations that links customer 
service level to arbitrary probability distribution functions (including empirical distributions) 
reflecting repair lead-time and part arrival at a repair shop. The method represents a customer 
service level as a mean of the number of on-time delivered parts during a unit time interval 
divided by a mean number of arrivals per time unit to obtain a distribution of on-time 
delivered parts per time unit. A search procedure is used to locate the optimal inventory level 
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using the distribution. Assets having multiple rotable part types can be analyzed using a 
constrained optimization method. 

[15] Figure 1 is a block diagram illustrating the rotable inventory problem to be addressed 

and optimized. As shown in the figure, rotable part arrivals 100 reach a repair shop according 
to an arrival probability distribution, which reflect the pattern in which parts arrive at a repair 
shop at a given time, and go through the repair process having a repair lead-time 102 with its 
own probability distribution. The arrival distribution and the repair lead-time distribution may 
be different from each other and, due to the nature of real world constraints, is uncertain. 

[16] The customer service level is determined based on a desired turn around time (TAT) 

104. The customer service level is defined as the probability that a given part repair 
turnaround time is less than or equal to the customer's desired/contracted turnaround time 104. 
In Figure 1, the path containing the desired turnaround time 104 corresponds to customer 
requests for obtaining repaired parts or replacements back. The desired turnaround time 
distribution 104 and the repair lead-time distribution 102 are then considered together to 
determine the optimum rotable inventory level 106 from which deliveries of repaired and 
replacement parts will be drawn. 

[17] The algorithms through which the optimal rotable inventory is calculated will now be 

described. For purposes of explanation in the equations below, I represents the rotable 
inventory level 106, W the desired turnaround time 104, and a represents a desired on-time 

delivery rate. The inventive method aims to find a minimum rotable inventory level I* at 
which the probability that customer turnaround time is less than or equal to the desired 
turnaround time value W is greater than or equal to the given desired probability level a . In 
other words, the method aims to find a minimum rotable inventory level I* at which there is a 
high probability (as defined by ot) that the inventory level will be high enough to meet or 
exceed customer service expectations with respect to the turnaround time. 
[18] As noted above, the repair lead-time is uncertain, which therefore makes the 

customer's turnaround time uncertain as well. In other words, there is no way to predict how 
quickly a customer can get a repair done or a replacement part installed With the desired 
turnaround time 104, customer service level is defined as the probability that the customer 
turnaround time for getting a repair or replacement is less than or equal to the 
desired/contracted turnaround time 104, that is, the probability that a customer part demand is 
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satisfied in a timely fashion. In the equation below, wtat represents the random customer 
turnaround time (i.e., the time needed for a given part repair or replacement). Let a represent 
the desired service level. Then the customer service level requirement can be expressed by the 
following equation: 



According to the mathematical Law of Large Numbers, which is known in the art, the 
probability described in Equation 1 can be rewritten as: 



where m# is the total number of on-time-delivered parts by time k (e.g., a day) and the total 

number of parts that need to be repaired (i.e., total repair demand) by time k. Equation 2 
holds true because the ratio of on-time delivered parts to the total number of repaired parts 
will reflect the degree to which customer-requested turnaround times have been met. 

If the number of parts arriving at the repair shop during a given time unit k is constant, 
the total number of parts that need to be repaired n* during the time unit can be expressed as: 



where A is the number of parts that arrive by time k. As can be seen in Equation 3, the total 
repair demand n k therefore is the number of parts A that arrive by time k multiplied by a 
difference between the time k and the desired turnaround time W. 

For the discussion below, assume that a random repair lead-time r is within a 
minimum repair time r and a maximum repair time R (i.e., in the interval [r, /?]). Based on 
this information, it is known with certainty that parts that are starting to be repaired at and 
before time (k - R) will be completed at time fc, and parts starting after time (k - r) will still be 
in repair at time k. Parts having repair start times that start within time [k - /?, k - r], however, 



Pv{w TAT <W)>a 



(1) 




(2) 



n k =A-(k-W) 



(3) 



EXPRESS MAIL NO. EV220475805US 67,097-043 

EH- 11224 

may or may not be done by time k; in other words, there is a probability, rather than a 
certainty, that these parts will be done by time k. The total number of repaired parts that will 
be completely repaired by time k can therefore be expressed as: 



R k =A(k-R) + X j^Xyik) (4) 

i=k-R+\ y=l 



[22] 



where X(j{k) is a random binary variable equal to either 0 or 1. Xyifc) represents a repair status 

at time k for a given f h part that arrived at time i, where j =7, 2, ...A; and i = k - R + 1, k- R 
+ 2, ...n - r. X(j(k) will equal 1 if the part is done being repaired at time k and will equal 0 

otherwise. 

Because A is a constant in this example, reflecting a constant number of part arrivals 
per time unit k, the number of on-time deliveries at a time k (i.e., OT*) will fall in the interval 
[0, A]. 



OT k = Min (A, Max (0, 1 + R k - n k _ x )) 



- Min 



A, Max 



k-r A 



0,/-A.(/?^W) + A+ X 

i=k-R+\ 7=1 J J 



(5) 



[23] In other words, the number of on-time deliveries will be the smaller of: (1) the total 

number of parts that have arrived A or (2) the larger of 0 and the difference between the total 
rotable inventory level and the number of parts whose repairs have been completed by time k. 

[24] The total number of on-time deliveries by time k can therefore be expressed as: 

1=1 



By replacing m£ and in Equation (2) with the expressions in Equations (3) and (6), the 
probability of on-time delivery of repaired parts can be expressed as: 
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Pr(w TAT <W)=lim 



— • lim 



A, 



or 



(7) 



where // or represents the mean on-time delivery of replacement/repaired parts. With the 
expression in Equation 7, it is now possible to find the minimum inventory level / (which 

determines the mean on-time delivery jx 0T ) such that > a (i.e., the mean number of on- 

A 

time delivered parts divided by the number of requested parts per time unit, which represents 
the actual delivery rate, is greater than or equal to the desired on-time delivery rate a). 
Expressed another way, it is now possible to find the minimum inventory level I* such that 
Mot Za-A. 

[25] For this calculation, assume that Xy(k) for all (i, j) are independent; that is, the repair 

of one part is not affected by the repair of another part. According to the Central Limit 

k-r A 

Theorem, which is known in the art, the term z k = Z ^X 0 (&) used in calculating the 

number of on-time deliveries at time k (OT*) in Equation 5 above has a normal distribution 
with the following mean jll and variance a. 

V = Mk = £ Z E ( X iM))= A ' Z A- £ G, (8) 

i-k-R+\ j-l i=k-R+\ i=r 

° 2 =°l= Z G ( (l-G,) (9) 

i = jt-/? + l 7 = 1 f'= r 

where G; represents the repair lead-time distribution (i.e., the probability that repair lead-time 

will be less than or equal to time i). 
[26] Referring to Figure 2, it can be seen that the term (/ + /?£ - n^j), which reflects a 

given inventory level (I) plus the total number of parts repaired by time k (i.e., R*) minus the 
total demand at time M (i.e., n*.y), has a normal distribution O , as represented by Series 2 in 
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Figure 2, with a mean of = I -A (R-W -/)+// and a variance a 2 . As also shown in 
Figure 2, the on-time delivery OT^ has a distribution *P , as represented by Series 1 in Figure 

2, that is truncated from the normal distribution O by removing the tails that lie beyond the 
area bounded by [0, A] and adding the mass probability in each tail to the mass probability at 
the corresponding boundary for each tail. 
[27] From this information, it is simple to compute the mean on-time delivery ju 0T , and 

therefore the service level Pr {w TAT <w), that is possible for a given inventory level I. Given 
an inventory level I, the mean and variance of the normal distribution <D can be computed 
using Equations (8) and (9). The probability mass distribution ¥ of OTfc can also be easily 

computed by adding the probability in the tail portion of its distribution to the probability at 
the boundary corresponding to the tail. With this distribution Y , the mean on-time delivery 
ju 0T can be computed. 

[28] Computing the required inventory /* from a specified probability for satisfying a part 

demand within a given turnaround time (i.e., Pr(w TAT <w)), however, is less straightforward 
than computing the mean on-time deli very ju 0T because the non-linear minimum/maximum 
function used to define on-time delivery, as shown in Equation 5, has an irregular distribution. 
Because of this, a search procedure, such as the one shown in Figure 3, may be used to 
determine a required minimum inventory level I* that can fulfill a given desired turnaround 
time with a selected probability. 

[29] The method shown in Figure 3 assumes that the part repair lead-time distribution Gi is 

given. From this information, the method first computes the mean ju and variance a of the 
part repair lead-time distribution using Equations (8) and (9) (block 150). 

[30] Next, the method defines a probability term /? = Pr(g > A) where g has a normal 

distribution 3> as represented in Series 2 in Figure 2 (block 152). In other words, (3 reflects 
the probability that g is greater than or equal to the number of parts arriving per unit time. 
Because p is expressed as a probability, it is clear that /?e (o,l). Further, as reflected in the 
expression ju^ = I -A-(R-W -l)+ju noted above, //^ is a function of the inventory level /. 
The method initializes (3 as a (i.e., the desired on-time delivery rate) and sets its initial left 
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and right bounds as BL and BR, respectively (block 154). The values for BL and BR. can 
cover any desired range, and in one embodiment (3 lies midway between BL and BR. 
[31] Next, a difference term A^ is searched or computed from a standard normal 

distribution table reflecting the distribution of (5 via any known means (block 156); this 
difference term can be obtained by, for example, a standard normal distribution table. The 
inventory value I corresponding to the distribution of |3 can be calculated 
from/ = A'(R-W)- fi + A^ <r, and this corresponding inventory value I is treated as a 

proposed inventory value for computing the on-time delivery probabilities that are eventually 
used to generate a mean on-time delivery value to be compared with the desired mean on-time 
delivery value. 

[32] From the known probability mass distribution V of OTfc and the proposed inventory 

level I corresponding to |3, the method computes on-time delivery probabilities for OT =0, 1, 
2... A and determines the mean on-time delivery ju 0T based on the resulting on-time delivery 
probability distribution for each value of OT (block 158). 
[33] The mean on-time delivery value mot ls then compared with a desired mean on-time 

delivery value (i.e., the method checks if ju 0T <a-A) (block 160). If the mean on-time 
delivery value is less than the desired mean on-time delivery value, it indicates that the value 
selected for (3 and its corresponding distribution are too far to the left on a distribution graph, 
causing the corresponding inventory level I to be too low with respect to the desired 
turnaround time (e.g., the proposed inventory level I is insufficient). In response, the left 
bound BL is set equal to the proposed value (3 ar\dj3 = ± (BL+BR) (block 162). In other words, 

(3 and BL are reset to cut the distance between the left and right bounds in half and move the 
position of 3 to a new location midway between the two bounds. The process then returns to 
block 156 to recompute the probability distribution of f} with the new value for (3 as well as 

the new left and right bound values BL and BR. 
[34] If, however, the mean on-time delivery value is greater than or equal to the desired 

mean on-time delivery value, the difference between the mean on-time delivery value and the 
desired mean on-time delivery value is then compared with a selected convergence threshold 
e (i.e.,// OT -a-A<€ ) (block 164). The value selected for e can be any desired small value 
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reflecting convergence between the mean on-time delivery value and the desired mean on- 
time delivery value. If the difference between the mean and desired mean on-time delivery 
values is less than the convergence value, it indicates that the proposed inventory I calculated 
from the value for P and its corresponding value for Ap will result in a mean on-time delivery 
value that is nearly equal to the desired mean on-time delivery value within an acceptable 
degree. At this point, the calculated proposed inventory level reflects the minimum inventory 
level I* needed to fulfill the desired on-time delivery, and as noted above can be computed as 
/ = A • (R - W) - fi + Ap • a (block 166). 

[35] If the difference between the mean and desired mean on-time delivery values is greater 

than the convergence value e, it indicates that the current value selected for p and its 
corresponding distribution are too far to the right on a distribution graph, causing the 
corresponding proposed inventory level I to be higher than needed to meet the desired 
turnaround time. In response, the method lets BR = fi, p = \\BL+BR) (block 168) before 

returning to block 156 to recompute the probability distribution of /? with the new value for 
P as well as the new left and right bound values BL and BR. The process continues to iterate 
until the value for P results in an inventory level that causes the mean on-time delivery value 
to be close to the desired mean on-time delivery value within a desired convergence 
z(^ 0T -aA<e). 

[36] As noted above, the process shown in Figure 3 assumes that parts arrive at a repair 

shop in a constant manner, where the number of arrived parts A remains the same for each 
time unit (e.g., per day). The theory and calculations in a case assuming that part arrivals are 
random are somewhat different. Generally, the inventive process handles random part arrivals 
by finding a way to approximate the random arrivals with a constant arrival expression, 
thereby allowing the process shown in Figure 3 to be carried out for random part arrivals. 

[37] More particularly, if it is assumed that the number of parts arriving at a repair shop at 

time k is random and represented by a random variable A^, the total part demand by time k 
can be expressed as: 
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[38] The total number of repaired parts by time k, denoted as /fy> is also a random variable 

and includes repaired parts that are completed with certainty as well as parts that have a 
probability of being complete at time k. As explained above, given that R is a maximum 
repair lead-time and r is a minimum repair lead-time, it is known with certainty parts that 
arrive at and before time (k - R) are completed by time k, and parts starting after (k - r) are 
still in repair by time k. Parts that arrive within [k - /?, k - r] may or may not be done by time 
fc, causing them to be completed by time k with probability rather than with certainty. The 
total number of completed parts by time k can therefore be written as 

k-R R-\ 

** = Za+2Xa) (id 

t=0 t=r 



Similar to the example shown in Equation 5, the number of on-time deliveries at k can 
therefore be expressed as 



OT k = Max (0, Min {A k _ w , / + R k -n k _ x )) 



= Max 



-Max 



0, Min 



k~W~\ R-\ 



t=k-R+\ t=r 



0, Min 



k-W-\ 



J J 

k-r 



(12) 



4k-w.'+ 5>,(G ik _,-l)+ 5>,G*_, 

t=k-R+\ t=k-W 



J) 



The problem is now to find the mean on-time delivery value Hot that will be greater than or 
equal to a desired mean on-time delivery value (i.e., ju 0T >a A), which will cause the 
inventory level calculated from the mean on-time delivery rate \Xqt to be the minimum 
inventory level /* 

[39] Assuming that for all k are independent and with reference to the examples in 

Equations 8 and 9, the term z* used in calculating the number of on-time deliveries at time k 
can be expressed according to the Central Limit Theorem by 

k-W~\ k-r 

z k = £/*,(G(fc-/)-l)+ ^A t G(k-t) and will have a mean and variance of: 

t=k-R+\ t=k-W 
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k-W-\ k-r 



t=k-R+\ t=k-W 
( R-\ W \ 



(13) 



Jfe-W-l k-r 



<r 2 = °k= £M4C*J 



f=*-fl+l t=k-W 
(W R-\ \ 



(14) 



As noted above with respect to the constant daily arrival example, (/ + - n^j) also 

has a normal distribution 4> (Series 2 in Figure 2) with a mean of and a variance a 2 . 
Because A# is random in this example, the distribution of OT^ is complicated. The 
distribution of OT^ can be numerically calculated. The mean fi 0T can then be computed. 

More particularly, assume that an arrival value A^y/ and an inventory value 

k-W-\ k-r 

/+ ^ A t( G k-t ~1)+ ^ A t G k-t are independent with their own respective distributions G and 

t=k-R+\ t=k-W 

O , respectively. The distribution of the Max ( • ) term in OTfc can then be numerically 
calculated as: ). 



By forcing w > 0 in the above distribution and adding the mass probability in the tail to the 
probability at 0, the distribution of OT is obtained in the same manner as in the constant 
arrival case. 

If the variance in the random part arrivals is low, it is possible to approximate the 

random arrival with a constant daily arrival term A . The algorithm for constant arrival can 
then be used even where the arrivals are actually random. 



^max (w) 55 G(w) + - G(w)- O(w) 



(15) 



12 



EXPRESS MAIL NO. EV220475805US 67,097-043 

EH- 11224 

[43] The examples above focus on optimizing inventory for a single part, but the inventive 

method can be applied for assets having multiple rotable parts without departing from the 
scope of the invention. Existing methods for constrained optimization problems can be 
applied to the inventive method to handle interrelations among rotable parts to find the 
optimal inventory for each part in the asset. 

[44] The algorithm described above can be extended to calculate optimal rotable inventory 

levels for individual parts in a single asset. Given an expected service level for an asset, the 
algorithm can determine optimal service levels and inventory levels for individual parts to 
fulfill the expected service level for the entire asset. In this example, it is assumed that an 
asset has a total of / rotable parts, indexed as i =1, 2, The expected service level for the 
entire asset is p. With the service levels for individual parts denoted as p\,P2>—Pi > the 
following condition needs to be satisfied to meet the expected service level for the entire 
asset: 

/ 

UPi^P (16) 
i=l 

[45] The left side of Equation 16 is the asset service level that is determined by the service 

levels of individual parts. As shown in previous examples, given a service level p t for a part 

i, the minimum rotable inventory can be computed for that part. 
[46] More particularly, let jc,- represent the computed inventory level, and let SI represent 

the relationship defined by the algorithm illustrated in previous examples. The computed 

inventory level jt f - can be expressed as: 

x t = MPi) (17) 

Let c x be the cost for having one unit of part i in inventory. Thus, the total rotable inventory 
cost can be written as: 

/ 

J = "£ x i' c i (18) 
*=1 
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The rotable inventory levels for individual part types within the asset can therefore be 
obtained by solving the constrained optimization problem defined by Equations 16, 17 and 18. 
That is: 

/ 

min / with / = £jc £ q (19) 

subject to 
/ 

Upi^p 

1=1 

*i = MPi)> i=l,2, ; ../. 
ft 6 (0,1], i=l,2,.../. 

[47] With this constrained optimization problem defined, existing constrained optimization 

methods can be used to obtain the optimal rotable inventory levels for each part type in an 
asset to meet the customer service level for that asset. 

[48] The minimum inventory value I* therefore reflects the optimum rotable inventory 

level for a given service level, a given repair arrival, and repair time statistic data. By 
generating a service level and checking whether the average service level meets a desired 
service level, it is possible to calculate the minimum rotable inventory level that can satisfy 
the desired service level. 

[49] Figure 4 is a block diagram of a computer system 200 that can be used to implement 

the method according to one embodiment of the invention. The system includes a user 
interface 202, a processor 204 that executes software that can carry out the inventive 
algorithm, and a memory 206 that stores data such as look-up tables, part type characteristic 
data, turnaround time information entered by the user, etc. The user interface 202 can be any 
system that allows a user to send data to and receive data from the processor 204, such as a 
keyboard/monitor combination, touch screen, graphical user interface, etc. 

[50] It should be understood that various alternatives to the embodiments of the invention 

described herein may be employed in practicing the invention. It is intended that the 
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following claims define the scope of the invention and that the method and apparatus within 
the scope of these claims and their equivalents be covered thereby. 
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